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Abstract 

The ground state hyperfine splitting values of high Z hydrogenlike 
ions are calculated. The relativistic, nuclear and QED corrections are 
taken into account. The nuclear magnetization distribution correc- 
tion (the Bohr-Weisskopf effect) is evaluated within the single particle 
model with the gg-factor chosen to yield the observed nuclear moment. 
An additional contribution caused by the nuclear spin-orbit interac- 
tion is included in the calculation of the Bohr-Weisskopf effect. It is 
found that the theoretical value of the wavelength of the transition 
between the hyperfine splitting components in 165 ^ 66 + i s i n good 
agreement with experiment. 

PACS number(s): 31.30.Gs, 31.30.Jv 
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1 Introduction 



Laser spectroscopic measurement of the ground state hyperfme splitting in 
hydrogenlike 209 Bi 82+ [1] has triggered great interest to calculations of this 
value (see [2-7] and references therein). Recently the transition between the 
F = 4 and F = 3 hyperfme splitting levels of the ground state in hydrogenlike 
i65j^ 66+ wag ODservec i [3] anc i jt s wavelength was determined to be 572.79(15) 
nm. It was found [8] that this value is in disagreement with commonly 
tabulated values of the nuclear dipole magnetic moment of 165 Ho (see, e.g., 
[9]) and in good agreement with the value that was measured by Nachtsheim 
[10] and compiled by Peker [11]. 

The hyperfme splitting values of 165 Ho m+ (for the magnetic moment from 
[9]) and other possible candidates for such experiments were evaluated, with- 
out QED corrections, in [4]. In the present paper we refine the calculations 
of [4] considering a more accurate treatment of the nuclear effects and taking 
into account the QED corrections. 



2 Basic formulas and calculations 

The ground state hyperfme splitting of hydrogenlike ions is conveniently writ- 
ten in the form [4,12] 



A£ M = -a(aZ) 



4 3 ji m 21 + 1 



rnc 
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3 /in m>p 2/ 

x{A(aZ)(l-5)(l-e)+x rad } . (1) 

Here a is the fine structure constant, Z is the nuclear charge, m is the electron 
mass, m p is the proton mass, \x is the nuclear magnetic moment, fi^ is the 
nuclear magneton, I is the nuclear spin. A(aZ) denotes the relativistic factor 
[13] 

• 4 < aZ ) = ^TTT) = 1 + f (aZ > 2 + T (aZf + " ' • (2) 



where 7 = yl — (aZ) 2 . 5 is the nuclear charge distribution correction, e is 
the nuclear magnetization distribution correction (the Bohr-Weisskopf cor- 
rection) [14], and x rad is the QED correction. 
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To calculate the nuclear charge distribution correction 5 we used the two- 
parameter Fermi model: 



p( r ) = £ (3) 

I y ! 1 + exp [(r - c)/a] V ; 

The parameters c, a and (r 2 ) 1 / 2 have been taken from [15]. We found that 
the error due to an uncertainty of the nuclear charge distribution parameters 
is much smaller than an uncertainty of the Bohr-Weisskopf effect. The Bohr- 
Weisskopf effect was calculated assuming that the magnetizations can be 
ascribed to the single particle structure of the nucleus with the effective gs- 
factor chosen to yield the observed nuclear magnetic moment. The nucleon 
wave functions were calculated by the Schrodinger equation with the Woods- 
Saxon potential [16-17] 

U(r) = V(r) + V so (r) + V Cou i(r) , (4) 

where 

V(r) = -V f(r) , 
V so = \Vv{h/2m p c) 2 o- ■ 1 r- l df so {r)/dr , 
a(Z- l)(3-r 2 /i? 2 )/2i? , r < R \ 



Vcoul — 



a(Z-l)/r, r>R 



f(r)= l + exp[(r-R )/a] 



-l 



fso(r) = [1 + exp [(r - R S o)/a] 

In the neutron case, the term Vcoui must be omitted. The phenomenolog- 
ical spin-orbit (SO) interaction, characterized by the parameter A, is much 
larger than it follows from the Dirac equation. In the present paper we 
used the potential parameters from [17]: Ro = 1.347-As fm, Vq = 40.6 MeV, 
R so = 1.280A5 fm, A = 31.5 for neutron and R = 1.275A5 fm, V = 58.7 
MeV, Rso = 0.932v43 fm, A = 17.8 for proton. Despite these parameters 
were chosen to give reasonable binding energies for the lead region alone, we 
used them at lower Z as well. We found that the uncertainty of the Bohr- 
Weisskopf effect caused by possible changes of these parameters [18,19] is 
small in comparison with an expected error due to deviation from the single 
particle model. 
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As is known [20], the SO interaction gives an additional contribution 
to the nuclear magnetic moment. Taking into account the related term in 
the hyperfine splitting theory gives an additional contribution to the Bohr- 
Weisskopf effect. If we denote the SO interaction by 



V SO {r)=<i>so(r)(s-l), 



(5) 



the total expressions for the Bohr-Weisskopf correction within the single par- 
ticle approximation are given by 



9s 



9i 



1 



(27 - 1) 



+ 



9l_ 
9i 



(21 



21 



■(K L ) + 



(21 + 1) m p 
47(7 + 1) h 2 



(<Psor 2 K L ) 



for 7 = L + \ and 



9s_ 
9i 



2(7 



(Ks) 



(27 + 3) 
87(7+1) 



(K s - K L ) 



+ 



9l 



9i 



(27 + 3) (27 + 1) m Vlx 2 

L ' ~ TrrrrrtTTWsor k l ) 



for 7 



Here 



2(7 + 1) 

(Ks) 
(Kl) 



(<Psor 2 K L ) 



47(7 + 1) h 1 



Ks(r)\u(r)\ 2 r 2 dr , 

KL(r)\u(r)\ 2 r 2 dr , 
<P so (r)r 2 K L (r)\u(r)\ 2 r 2 dr , 



K s (r) 



Jo fg dr' 
fo°°f9dr> ' 



S r Ai- r 4)fgdr> 



(6) 



(7) 



Kl{r) JTfgdr' ' 

g and / are the radial parts of the Dirac wave function of the electron, u(r) 
is the radial part of the wave function of the odd nucleon. The functions 



4 



Ks(r) and Kl(t) are calculated by using simple approximate formulas from 
[4] (the relative precision of such a calculation is of order aZR /(h/mc)). 
Taking g L = for neutron and Ql = 1 for proton, we choose gs to give 
the experimental value of the magnetic moment within the single particle 
approximation: 

A 4 1 \t 1 27 + 1 m p . I 2 ,-| /o , 
for 7 = L + | and 

7(27 + 3) 27+1 m p ,, 



li N 2(7+ L 2(7 + 1) 4(7+1) ft 2 



(</>soO <?l (9) 



for 7 = L - i. 

In the third and fourth columns of the table 1 we present the values 
gs and e calculated by equations (6)- (9). As one can see from the table, 
except 159 T6 and 127 7, the values gs lie between the free Dirac and free real 
g-factors. For comparison, in the fifth column we give the values e found 
in disregarding the spin-orbit terms in equations (6)-(9) (it corresponds to 
the calculation using the original Bohr-Weisskopf formulas [14]). In the last 
column we give the values e found in [4] by using a simple, homogeneous 
over the nucleus, distribution of |w(r)| 2 . (We note here that in the case of 
209 Bi S2+ in [4] a more accurate evaluation of e was also presented which 
gave e = 0.013 and A = 242.0 nm, without the QED correction. The same 
value was found in [21].) Taking into account that the single particle model 
with the effective gs-factor gives reasonable agreement with experiment for 
neutral atoms [22,23], we assume the following errors bars for e. For the ions 
where the spin and orbital parts in (6)- (7) are of the same sign (it results 
in relatively large values of e) the uncertainty is about 30% of e. For the 
ions where the spin and orbital parts in (6)- (7) are of the opposite sign (it 
results in relatively small values of e) the uncertainty is about 20% of (Kg) ( 
0.2(7\~s) ~ 0.03aZb, where b is a factor tabulated in [4]). In the case of Pb the 
uncertainty is assumed to be about 10% of e. It should be stressed however 
that the uncertainty found in this way is to be considered only as the order of 
the expected error. More accurate calculations of the Bohr-Weisskopf effect 
must be based on the many particle nuclear models and must include a more 
consequent procedure for determination of the error bars. 
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The radiative correction is the sum of the vacuum polarization (VP) and 
self energy (SE) contributions. The VP contribution can easily be calculated 
within the Uehling approximation. We calculated this effect for a finite 
nucleus charge distribution and found that for Z=82, 83 our results are in 
good agreement with the results of [3,6]. The values of the electric loop 
and magnetic loop contributions in the Uehling approximation are given in 
the second and third columns of the table 2. The Wichman-Kroll (WK) 
contribution is also the sum of two terms. The first term is given by the 
WK electric loop correction to the electron wave function. The calculation 
of this term can be done in the same way as the calculation of the first 
order WK contribution. The results of such a calculation, based on using the 
approximate formulas for the WK potential for a point nucleus [24], are given 
in the fourth column of the table. The second term is the WK magnetic loop 
contribution. Calculation of this term is a more complex problem. However, 
calculations of the corresponding term in the VP screening diagrams for 
two-electron ions [25] allow us to expect that this term is small enough. In 
the fifth column of the table 2 the total VP contribution, without the WK 
magnetic loop term, is given. The SE contribution was evaluated in [6,26] 
in a wide interval of Z for a finite nuclear size distribution. The values 
of this contribution given in the sixth column of the table 2 are found by 
interpolation of the related values from [26]. In the last column of the table 
the total radiative correction is listed. 

In the table 3 we give the theoretical values of the energies and the 
wavelengths of the transition between the ground state hyperfine splitting 
components of high Z hydrogenlike ions. The error bars given in the table 
are mainly defined by the uncetainty of the Bohr-Wesskopf effect discussed 
above. The magnetic moment values, except Ho [10,11,8], are taken from[9]. 
In the case of Ho, the theoretical value constitutes A = 572.5(1.7) nm and 
is in good agreement with experiment A = 572.79(15) nm [8]. The values 
of the individual contributions are listed in the table 4. In the case of Bi 
the difference between the experimental value (A = 243.87(2) nm [1]) and 
the theoretical value given in the table 3 is within the expected error of the 
Bohr-Weisskopf effect. We expect that these theoretical results will be re- 
fined by including the RPA approximation in a more elaborate treatment of 
the Bohr-Weisskopf effect [5], based on the dynamic-correlation model [27]. 
Such calculations are underway and will be published elsewhere. 
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Table 1: The Bohr-Weisskopf effect within the single particle model of the 
nucleus: with taking into account the SO term in (6)- (9), without the SO 
term in (6)-(9), and taken from [4]. 



ion 


Nucleon state 


9s 


e (with SO) 


e (without SO) 


e (from [4]) 


US /n 48+ 


Iff* 


3.67 


0.0047 


0.0045 


0.0039 


im 56 50+ 


2d 5 


3.04 


0.0052 


0.0051 


0.0046 


1113 S6 60+ 


l n 


4.21 


0.0014 


0.0019 


0.0016 




2d 5 


1.95 


0.0052 


0.0051 


0.0047 


WA Cs U+ 




4.14 


0.0017 


0.0024 


0.0020 




Iff? 


3.64 


0.0025 


0.0032 


0.0027 


141p r 58+ 


2d 5 


4.88 


0.0075 


0.0073 


0.0072 


151£ U «U+ 


2d5 


3.27 


0.0080 


0.0079 


0.0079 


159^^64+ 


2d 3 


-0.203 


0.0069 


0.0074 


0.0073 


165^J 6b+ 




2.90 


0.0089 


0.0085 


0.0086 


175 Lu 70+ 


Iff? 


5.10 


0.0006 


0.0020 


0.0018 


l8l Ta 7^+ 


l5i 


4.76 


0.0017 


0.0032 


0.0030 


185^74+ 


2d 5 


2.71 


0.0122 


0.0120 


0.013 




3s 1 


3.47 


0.0179 


0.0177 


0.020 


2U5j>;8U-f- 


3si 


3.50 


0.0179 


0.0177 


0.020 


2()7p^81 + 


3pi 


-3.56 




0.0419 


0.036 




1/19 

■* 


2.80 


0.0118 


0.0133 


0.011 
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Table 2: The radiative corrections to the ground state hyperfine splitting 
in terms of x defined by equation (1). Xyp IS the Uehling electric loop 
contribution, X y p IS the Uehling magnetic loop contribution, x v / p E is the 
WK electric loop contribution, xyp is the total VP contribution without the 
WK magnetic loop part, xse is the SE contribution found by interpolation 
of the related values from [26], and x ra d is the total radiative correction. 



z 


~UE 

Xyp 


Xyp 


Xyp 


xvp 


XSE 


Xrad 


49 


0.0020 


0.0011 


-0.0000 


0.0031 


-0.0074 


-0.0043 


53 


0.0024 


0.0013 


-0.0000 


0.0036 


-0.0084 


-0.0048 


57 


0.0029 


0.0014 


-0.0000 


0.0043 


-0.0096 


-0.0053 


63 


0.0037 


0.0017 


-0.0001 


0.0054 


-0.0116 


-0.0062 


67 


0.0044 


0.0020 


-0.0001 


0.0063 


-0.0132 


-0.0069 


71 


0.0053 


0.0022 


-0.0001 


0.0074 


-0.0150 


-0.0076 


75 


0.0064 


0.0026 


-0.0002 


0.0088 


-0.0171 


-0.0083 


82 


0.0089 


0.0033 


-0.0003 


0.0119 


-0.0218 


-0.0099 


83 


0.0094 


0.0034 


-0.0003 


0.0125 


-0.0226 


-0.0101 



Table 3: The energies (AE) and the wavelengths (A) of the transition be- 
tween the hyperfine structure components of the ground state of the hydro- 
genlike ions. A is the relativistic factor, 5 is the nuclear charge distribution 
correction, e is the nuclear magnetization distribution correction (the Bohr- 
Weisskopf correction), and x ra d is the radiative correction see equation (1)). 



Ion 


M 


A 


6 


e 


^rad 


AE (eV) 


A (nm) 


113y- n 48+ 


5.5289(2) 


1.2340 


0.0170 


0.0047 


-0.0043 


0.9148(13) 


1355.3(1.9) 


V2L sb M+ 


3.3634(3) 


1.2582 


0.0191 


0.0052 


-0.0045 


0.6891(11) 


1799.3(2.8) 


L'2iS 


2.5498(2) 


1.2582 


0.0191 


0.0014 


-0.0045 


0.4994(7) 


2482.5(3.5) 


L'27 I b'2+ 


2.81327(8) 


1.2843 


0.0213 


0.0052 


-0.0048 


0.6587(10) 


1882.2(3.0) 




2.58202 


1.3125 


0.0237 


0.0017 


-0.0051 


0.6582(11) 


1883.6(3.0) 


i3y La bb+ 


2.78305 


1.3430 


0.0263 


0.0025 


-0.0053 


0.8052(14) 


1539.9(2.6) 


141p r 58+ 


4.2754(5) 


1.3761 


0.0292 


0.0075 


-0.0056 


1.464(3) 


847.0(1.9) 


151 Eu 6'2+ 


3.4717(6) 


1.4509 


0.0365 


0.0080 


-0.0062 


1.513(4) 


819.4(2.0) 




2.014(4) 


1.4933 


0.0407 


0.0069 


-0.0065 


1.099(3) 


1128(3) 




4.132(5) 


1.5395 


0.0456 


0.0089 


-0.0069 


2.166(7) 


572.5(1.7) 


ivb Lu vu+ 


2.2327(11) 


1.6453 


0.0575 


0.0006 


-0.0076 


1.482(4) 


836.6(2.4) 


LiSl Ta T2+ 


2.3705(7) 


1.7061 


0.0645 


0.0017 


-0.0080 


1.758(5) 


705.3(2.2) 


lSb Re '/4+ 


3.1871(3) 


1.7731 


0.0706 


0.0122 


-0.0083 


2.749(10) 


451.0(1.7) 




1.62226 


2.0217 


0.0988 


0.0179 


-0.0096 


3.229(18) 


384.0(2.1) 




1.63821 


2.0217 


0.0989 


0.0179 


-0.0096 


3.261(18) 


380.2(2.1) 


M7 Pb ai + 


0.592583(9) 


2.0718 


0.1049 


0.0419 


-0.0099 


1.215(5) 


1020.5(4.5) 




4.1106(2) 


2.1250 


0.1111 


0.0118 


-0.0101 


5.101(27) 


243.0(1.3) 
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Table 4: The individual contributions to the ground state hyperfine splitting 
in 165 ff 66 + for 11 = 4.132(5)/Xjy [10,11,8]. 



Nonrelativistic value 


1.4945(18) eV 


Relativistic value (point nucleus) 


2.3007(28) eV 


Nuclear size effect 


-0.1050(7) eV 


Bohr-Weisskopf effect 


-0.0195(59) eV 


Vacuum polarization 


0.0094 eV 


Self energy 


-0.0197 eV 


Total theoretical value 


2.1659(66) eV [A = 572.5(1.7) nm ] 


Experiment [8] 


2.1645(6) eV [A = 572.79(15) nm ] 
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